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1 Introduction 

The evolution of the universe is well described by the standard hot big bang 
model and its early stage offers a common field for particle physics and 
astrophysics Not only have various observations in laboratories been clues 
to understand what occurred in the early universe, but new ideas have also 
led to resolutions of the problems concerning very early universe; e.g., grand 
unified models with strong first-order phase transitions led to the idea of 
inflation, and dark matter may be explained by new physics, which will be 
examined in the near future. On the other hand, astrophysical observations 
give constraints on models of particle physics which cannot be obtained from 
accelerator experiments 0]. The baryon asymmetry of the universe (BAU) is 
among these topics. 

The BAU is one of the most obvious facts|l|]: we do not observe any 
antimatter in our solar system, and high energy cosmic rays, which contain a 
small amount of anti-protons consistent as secondary products, are evidence 
of the baryon asymmetry on the galactic scale. The absence of hard 7 rays, 
which would be emitted on nucleon-antinucleon annihilation, from nearby 
clusters of galaxies, such as the Virgo cluster, implies that such a cluster 
consists of (1 ~ 100)Mgaiaxy — lO^"^""^"^ Mq of either baryons or antibaryons 
only. While we have no evidence for baryon asymmetry on a larger scale, 
matter should be separated from antimatter on scale of, at least, lO^^M©. 
We characterize this manifest quantity by the ratio of the baryon number to 
entropy 

ri£ ^ Hb - ni 

s s ' 

where s is the entropy density and nh{ni) is the (anti)baryon number density. 
This remains constant in the absence of a baryon-number-changing process 
and entropy production, during the expansion of the universe. To explain 
the light-element abundances within the framework of the standard big-bang 
nucleosynthesis, it is required that r] = ns/n^ = (1.5 — 6.3) x 10^^°p|. As 
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shown in Appendix A, the photon density is related to s by s = 7.04n^ 
at present]!]], so that 

— = (0.21-0.90) X 10^1°. (1.2) 

s 

This is very small but sufficient to form the present visible matter. Although 
one might think that matter and antimatter were separated on a scale of 
IO^^Mq in locally baryon-symmetric universe, the present BAU cannot be 
explained without any primordial asymmetry. The present baryons are those 
which survived the nucleon-antinucleon annihilation and which froze out at 
T ^ 20MeV. Since starting from the baryon-symmetric universe n^/s ~ 
ni/s ~ 7 X 10"^'^ after the freeze out, the BAU, which is about nine orders 
larger, could not be generated |jl]]. To avoid this annihilation, suppose that 
some hypothetical mechanism separated nucleons and antinucleons before 
T ~ 38MeV, when Ub/s = ni/s ~ 8 x lO"^"*^. At that time the matter 
contained in the causal region was only about IO^'^Mq <: IO^^Mq. Hence 
it is natural to assume that the universe had baryon asymmetry before it 
cooled down to T ~ 38MeV. 

To obtain the BAU starting from the symmetric universe, three condi- 
tions, ffist proposed by Sakharovf^, must be satisfied: 

(1) baryon number violation, 

(2) C and CP violation, 

(3) departure from equilibrium. 

Without condition (1), the symmetric universe remains baryon-symmetric, 
and it seems difficult to realize a locally asymmetric but globally symmet- 
ric universe, as we saw above. If C and CP were conserved, baryon num- 
ber could not be generated in the symmetric universe. This is understood 
as follows. Suppose that the initial state of the universe is described by 
a density operator po which is C- and CP- invariant. Since baryon num- 
ber is odd under C and CP, (^5)0 = Tr[po^B] = 0. If the hamiltonian 
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of the system is C- or CP-invariant, the density operator p at later time, 
whose time evolution is governed by the Liouville equation, is also invari- 
ant under C or CP transformation: [p, C] = or [p, CV] = 0. Because of 
CBC^^ = -B and [B,V] = 0, (ub) = Trfpn^] = TtIpCubC-^] = -Trfpn^], 
or {ub) = TT[pCVnB{CV)^^] = —TtIpub]. Here C(V) is the operator rep- 
resenting the charge conjugation (space inversion). In any case we have 
(^_b) = 0. Hence both C and CP must be violated to have nonzero {ub)- If 
baryon-number-changing processes are in chemical equilibrium, the chemical 
potential for the baryon number vanishes so that the equilibrium distribution 
of baryons coincides with that of antibaryons: ra^ = rag. This implies that 
even if the universe was baryon-asymmetric at some time, the BAU van- 
ishes after the universe experienced an equilibrium era when baryon-number 
changing processes were in effect. 

As the models of particle interactions satisfying conditions (1) and (2), we 
now have electroweak theories, grand unified theories (GUTs), supersymmet- 
ric extensions of them and others. Among these, baryogenesis within GUTs 
was first extensively studied[^ |1|. In the framework of GUT-baryogenesis, 
the main process with baryon number violation is the out-of-equilibrium de- 
cay of the heavy bosons X, which are the gauge or Higgs bosons of mass 
mx > lO^^GeV. Suppose that X decays into two channels qq {AB = 2/3) 
and ql {AB = —1/3) with branching ratio r and 1— r, respectively. C and CP 
are violated if r is not equal to the branching ratio f of the process X qq. 
Then the expectation value of the change in baryon number in the decay of 
X-X pairs is (AB) = |r - |(1 - r) - |f + |(1 - f) = r - f. If C or CP is 
conserved, r = f so that B is not generated.^] At T ~ mx, the decay rate of X 
bosons is roughly given by — amx, where a = (7^/(47r) with g being the 
coupling constant, (a ~ 1/40 for the gauge boson, a ~ 10~^ ^ for the Higgs 
boson.) The Hubble parameter is if ~ /mpi, where g^, ~ 10^~^ is 

^Here we do not see spin and momentum in the final states, so that P does not affect 

r. 
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the effective massless degrees of freedom. (See Appendix A.) At tempera- 
tures near mx, To — H, so tliat the chemical equilibrium between the decay 
and the production of X and X no longer holds. The Boltzmann equations 
are used to estimate the generated baryon number quantitatively[^]. 

Since the weak interaction incorporates C and CP violation, it and its 
extensions may be candidates to offer the BAU if they fulfill the other con- 
ditions (1) and (3). It is well known that in the standard model B + L is 
violated by the axial U{1) anomaly. This anomalous process can convert pri- 
mordial L into B. Hence if L is generated by some L- violating interaction at 
the intermediate scale between the electroweak and GUTs scales, B would be 
left after the electroweak phase transition. This might be another candidate 
for the baryogenesis. If the _B-violating process occurs out of equilibrium at 
the electroweak scale, the BAU might be generated within the framework of 
the electroweak theory. This possibility is rather attractive, since it depends 
only on physics which could be tested by near future experiments. In this ar- 
ticle, we review the attempt to realize this idea — electroweak baryogenesis, 
focusing on the relation between the BAU and CP violation in the Higgs sec- 
tor of the extensions of the standard model. For earlier works on this subject, 
see the review article, Ref. 0. This paper is organized as follows. In § 2, we 
present the basic idea of the electroweak baryogenesis. In the following three 
sections, we study how it is realized in some detail. We summarize recent 
attempts to determine CP violation at the electroweak phase transition in 
§ 6. The final section is devoted to concluding remarks. 

2 Overview of the electroweak baryogenesis 

The axial anomaly in the electroweak theory nonperturbatively violates the 
sum of baryon number and lepton number B + L and its probability is en- 
hanced at high temperatures in the symmetric phase of the gauge symmetry, 
while suppressed at low temperatures in the broken phase, as we shall see 
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in the next section. Near the electroweak phase transition (EWPT) tem- 
perature ~ lOOGeV, the universe was expanding too slowly to make the 
anomalous process out of equilibrium. If the EWPT is first order accompa- 
nying formation and growth of bubbles of the broken phase in the symmetric 
phase, the B + L-violating process will deviate from equilibrium near the 
bubble walls. C and CP violation of electroweak theory may distinguish 
the effects of the expanding bubble on the fermions and antifermions leading 
to the generation of B + L. If the anomalous process is fully suppressed in 
the broken phase, the BAU is left after the EWPT. We shall refer to this 
scenario as 'electroweak baryogenesis'. To investigate the possibility of this 
idea, we need a wide range of knowledge in theoretical physics; estimation 
of the nonperturbative B + L violating ratio, finite-temperature field the- 
ory to determine the static properties of the EWPT, djTiamics of the phase 
transition and model building of electroweak theory consistent with present 
experiments. 

Before surveying each subject, we briefly introduce other attempts to pro- 
duce the BAU by use of the anomalous B + L violation. Since B + L rather 
than B is violated, B + L is washed out if the anomalous process was in equi- 
librium. This implies that we need primordial B — L to have the present BAU 
if the anomalous process had been in equilibrium until it froze out0. (We 
shall discuss this issue in the next section.) Then the GUTs which conserve 
B — L, such as the SU{5) model, would be useless. As the origin of nonzero 
primordial B — L, one may consider a GUT which does not conserve this 
quantity. This assumes new physics at the GUT scale (~ lO^^GeV), which 
could be tested by the proton decay experiments. The L-violation at the 
intermediate scale between the electroweak and GUTs scales with the Majo- 
rana neutrino might seed the primordial -B — L0. For the L- violating process 
not to erase L completely, some conditions on the masses are imposed, which 
could be checked by solar neutrino and other experiments. Supersymmetric 
extensions of the standard model contain the scalar superpartners with the 
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same internal quantum number as the quarks and leptons. These scalar fields 
could have nonzero expectation values along the 'fiat directions', which are 
typical in such models, leading to B and/or L violation at high temperatures. 
This mechanism is investigated in Refs. and would be a probable origin if 
supersymmetry is discovered in the experiments. In addition to these, topo- 
logical defects, such as strings and domain walls, formed at the electroweak 
scale, could create the baryon number when they move or decay [|10|]. 

The present BAU might be generated by one or some combination of 
these mechanisms, including the electroweak baryogenesis. In any case, the 
electroweak baryogenesis would be the last chance to generate the BAU if 
the anomalous process froze out after the EWPT, so that it might affect 
the baryon number already created until that time. Once a model of the 
electroweak theory is specified, the nature of CP violation and the EWPT 
are, in principle, known, then the BAU generated can be evaluated. When 
quantitative study of the electroweak baryogenesis is developed, one can se- 
lect a model to explain the present BAU and such a model can be confirmed 
by experiments in the near future. This is why much attention is paid to 
this subject and related subjects such as the finite-temperature phase tran- 
sition of the gauge-Higgs system and the dynamics of the first-order phase 
transition. 



3 Sphaleron process 

3.1 Anomalous fermion number nonconservation 

In the standard model and its extensions with the same fermion-gauge inter- 
action, the current of B + L, suffers from the axial anomaly||ll||: 

d,fB+L = ^,[g'Tt{F,^Fn-9"B,,B^''], (3.1) 
d,fB-L = 0, (3.2) 
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where Nj is the number of the generations, g {g') and F^^, (B^^) are the 
gauge couphng and the field strength of the SU (2) l {U (l)y) gauge field v4^(s) 
{B^{x)), respectively, and the tilde denotes the Hodge dual, F^'^ = ^e^'-'^'^F, 
Integrating these equations over {t,x) E [U^tf] ® R-^, we obtain 



pa- 



B{tf) - B{U 



327r2 



d X 



g'TT{F^,Fn - g' B^^B 



Nf[Ncs{tf)-Ncs{ti)] 



(3.3) 



where Ncs is the Chern-Simons number defined, in the Weyl gauge {Aq = 0), 
by 



Ncsit) 



1 



327r2 



(fxe 



ijk 



g'Tr F,,A, 



;gAiAjAf 



g'^BijBk 



The Chern-Simons number is not gauge-invariant but its difference 
apparently gauge- invariant. For a pure-gauge configuration, Fij = B 



(3.4) 

:! is 

0, 



Ncs takes an integer value. That is, the classical vacua of the gauge sector 
are labeled by the integer Ncs, which is the winding number correspond- 
ing to the homotopy group 7!'3{SU{2)) ~ Z, since the U{1) contribution 
in (|3.4| ) always vanishes for a pure gauge. In the semi classical approxi- 
mation, the vacuum-to-vacuum amplitude with ANcs = 1 in the SU{2) 
gauge system is dominated by the instanton configuration | p. 2| | , which is a 
classical solution with finite euclidean action Sinst = ^n'^/g^, and is given 
by ~ exp(— 5'inst) |T3|- For a strong-coupling theory such as QCD, this am- 
plitude is so large that the quantum vacuum is a superposition of these 



degenerate classical vacua known as the 6'-vacuum|T^. As for the SU{2) 
gauge-Higgs system, there is no classical solution with a finite euclidean ac- 
tion, but the constrained instanton[]I5| or the valley instanton []TB| plays a 
similar role and yields the same suppression ~ exp(— Sinst) together with 
other factors. For the standard model, this suppression factor is so small, 
exp(— Sinst) = exp(— 27r sin^ 6'iy/aem) — e~^^^ = 10~^^, that the vacuum-to- 
vacuum transition with ANcs = 1 is not observable. That is, B + L can 
hardly change in the vacuum at zero temperature. 
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Figure 1: Schematical view of the energy as a functional of the field config- 
uration. The horizontal axis denotes the Chern-Simons number, which is a 
subspace of the configuration space. 



The nonconservation of the current (|3.1| ) indicates that B + L associated 
with the quarks and leptons actually changes when the background gauge 
fields have nontrivial topology. The Atiyah-Singer index theorem says that 
the integral / (i'^Tr[FF] — the Pontrjagin index — is related to the zero 
mode of the chiral fermionp!7|. This can be seen more directly if one considers 
the Dirac equation in this background of the gauge fields. Regarding the 
euclidean time as a parameter, the solution to the Dirac equation for the 
left-handed fermion in the instanton background shows that the level in the 
Dirac sea moves to a positive-energy level, as the parameter goes from —oo 
to oo|T3]. This is known as the level crossing phenomenon. 



3.2 Sphaleron transition 

Our main concern is whether the anomalous process is enhanced at high 
temperatures. This possibility was first suggested when the energy barrier 
between the neighboring classical vacua labeled by Ncs was found to be finite. 
The top of the barrier in Fig. |1| corresponds to the saddle point in the con- 
figuration space, which is a static saddle-point solution with finite energy of 
the SU{2) gauge-Higgs system p!9|1. Such a solution is called a "sphaleron" . In 
contrast to a topological soliton such as nonabelian monopoles and Skyrmions, 
it is unstable since the fluctuation spectrum around the solution contains one 
negative mode. It is a nontopological solution and has Ncs = 1/2 hi some 
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gauge, in which the vacua divided by it have Ncs = and 1. The negative 
mode represents the instabihty against decay to the neighboring vacua. Al- 
though a 'sphaleron' is not a topological soliton, its existence is related to 
the topology of configuration space of a field theory; some class of field the- 
ories having 'noncontractible loop' or its higher dimensional generalization 
may have such a saddle-point solution [|l9i Among these theories are the 



1 + 1-dimensional U{1) gauge-Higgs system [P2| and the 1 + 1-dimensional 
0(3) nonlinear sigma model [|23|, both of which also have instanton solutions. 
The static energy of the sphaleron in the system of SU{2) gauge fields and 
a Higgs doublet is given by 

i?sph = ^5 f 41 - lOTeV, (3.5) 
aw \g^J 

where A is the Higgs self coupling, aw = fl'^/(4vr) and 1.5 < i? < 2.7 for 
\/ E [0, oo)pD[. The sphaleron solution in the standard model with nonzero 
9w has been found and its energy is somewhat smaller [El . 

For temperatures below the barrier height, the transition rate between the 



classical vacua was estimated semiclassically by Affleckj^, whose classical 
statistical version had been formulated by Langer|26||. Affleck showed, in 
the WKB approximation of a quantum mechanical problem with metastable 
potential, that for lower temperature than (27r), where uj^ is the negative- 
mode frequency at the top of the barrier, the transition rate is given by 
r ~ 2 ImF, while for T > cj-/(27r). 



~— ImF. (3.6) 



vrT 

Here the free energy F = — TlogTr[e~^/^], expressed in the path-integral 



form, should be estimated around the dominant configurations, the bounce|27 



at low temperatures and the top of the barrier (sphaleron) at high temper- 
atures. In any case, the 'imaginary part' arises from the Gaussian integral 
around each configuration, which contains a contribution from one negative 
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mode. At T = 0, in the limit that the metastable vacuum becomes degen- 
erate with the stable one, the bounce action approaches twice the instan- 
ton action, so that T ~ exp(— 2S'inst), reproducing the instanton calculation. 
In field theories, the Gaussian integral contains contributions from positive 
modes as well as zero modes corresponding to global symmetries, such as 
translation and isospin, which are violated by the configuration. 

The first estimation of the transition rate in the SU{2) gauge-Higgs sys- 



tem was made by Arnold and McLerran||2^ and is given by 



rSUW„A/-„.|;^(^)'e--'- (3.7) 

where awiT) is the temperature-dependent fine structure constant, A/tr (A/'rot) 
is the zero-mode contribution from translation (rotation), which is about 26 
(5.3 X 10^) in the case A = g^, and k denotes the contributions from the other 
modes. Here ~ (1.8 ~ 6.6)m^ for 10~^ < X/g"^ < 10 and k is shown to 
be 0(1)111. 

All of these results are valid in the existence of the mass scale set by 
the vacuum expectation value of the Higgs, v. At T > Tc, where Tc is the 
EWPT temperature, v{T) = 0, so that the above results can no longer be 
applied. For such high temperatures, we have no theory but on the grounds 
of dimensional analysis [^, we have 

ri% ^ KiawT)\ (3.8) 

Although there is no sphaleron solution in the symmetric phase, we use the 
terminology 'sphaleron transition' for the anomalous process. The formula 
( |3.8| ) was checked by classical Monte Carlo simulations. The outline of the 
simulation is as follows: First, an ensemble of the coordinate and momentum 
of the classical field on a cubic lattice is generated by usual Metropolis or 
heat bath methods, according to the statistical weight exp{—H/T), where H 
is the classical hamiltonian. To avoid the Rayleigh- Jeans instability peculiar 
to thermodynamics of classical fields, the finite lattice spacing is adjusted at 
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each temperature. Picking up one of the configurations in this ensemble as 
the initial condition, the classical equations of motion are solved numerically. 
Thus {Nl;g{t)) is measured and is fitted to the expression for the random 
walk, {NQg{t)) = 2rVt as t ^ oo. This program was executed for a 3 + 1- 
dimensional SU (2) gauge- Higgs system in the symmetric phase and produced 
the result k > 0.4|3^, and for an SU{2) pure gauge system |3T[], which may 



be good approximation of the symmetric phase of the gauge-Higgs system, 
and produced the result n = 1.09 ± 0.04.0 

3.3 Washout of 5 + L 

As we noted in the previous section, an unsuppressed sphaleron transition is 
needed for the electroweak baryogenesis, but any generated baryon number 
may be washed out if the EWPT is second order or the sphaleron transition 
does not decouple after the EWPT0. For this decoupling to occur, the 
sphaleron rate should be smaller than the expansion rate of the universe. As 
we shall see in the next section, the EWPT took place at T = Tc ~ lOOGeV. 
At this temperature, the Hubble parameter is given by (see Appendix A) 

H{Tc) ~ 1.7y^^ - lO-^^GeV, (3.9) 
mpi 

where (7* ~ 100 is the effective massless degrees of freedom at this tempera- 
ture. At T > Tq, since the sphaleron transition rate per unit time, which is 
Pgpij in (|3.8|) multiplied by the particle density at that time, is much larger 



sph 

than the expansion rate, 

Tsph ^ f^atyT ~ 10-^ » if (Tc) , (3.10) 

the B + L-changing process is in equilibrium in the symmetric phase. As we 
show later, if the Higgs mass is larger than some value, even the sphaleron 



^Whether the system is in the symmetric or broken phase can be checked by mon- 
itoring the expectation value of a gauge-invariant operator such as ^^n)^{n) and 
^^n)Uf,{n)<^>{n + fi). 
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rate in the broken phase is larger than the Hubble parameter, so that the 
primordial B + L is washed out. 

The relic baryon number after this washout is estimated by use of the rela- 
tions among chemical potentials for the particles in the electroweak theory [P^. 
The particle number density for each degree of freedom at equilibrium is re- 
lated to its chemical potential by 



n+—n. 



1 1 

6" T 



for bosons 



for fermions, 
(3.11) 

where Uk = ^/k"^ + m? and we used m/T <: 1 and /i/T <; 1. Since the 
elementary processes in the standard model are in chemical equilibrium at 
high temperatures, there are relations among the chemical potentials. For 
example, the equilibrium of the gauge interaction imposes = /^d^ = 
fJ'iL — f^i = At_ + /io, where fiw is the chemical potential for W~ , if^dL) fo'^ 
the left-handed up-type (down-type) quarks, fiiL,R {l^i) for the i-th generation 
charged lepton (neutrino), and ji^ (fio) for the charged Higgs 0~ (the neutral 
Higgs 0°).^ Consider the electroweak theory with Nf generations and m 
Higgs doublets. When the sphaleron process such as |0) ^ uidLdLVL is in 
equilibrium, Nf{nu^ + 2/irf^) + f^i = 0. Various quantum number densities, 
such as B, L, the electric charge Q, and the isospin / are expressed in terms of 
these chemical potentials of the particles. By use of the relations among the 
particle chemical potentials, one can express the quantum number densities 
in terms of fiw, yUw^, /^o and = J2ifJ'i- 

B = Nf{fiuL + f^un + /irfi + /idfl) = 4:NffiuL + '^Nf^w, (3.12) 
L = Y^{^i + ^iiL + ^liu) = ^^i + 2Nf^w - Nf^iQ, (3.13) 

i 

2 1 



Since the strong interaction is in equilibrium, the chemical potential for each flavor 
is independent of color. Further, the mixing among the quark flavors is assumed to be in 
equilibrium. 
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- '^if^iL + fJ'iR-2- 2fiw - 2mfi_ 

i 

= 2Nf^iuL-2^i-{ANf + A + 2m)^iw + {ANf + 2m)^io, (3.14) 
11 1 

i 

= -i2Nf + m + 4)fiw, (3.15) 

where we normalized the densities by In the symmetric phase (T > 

Tc), all the gauge symmetries are respected so that Q = I3 = 0, which 
constrain two of the remaining chemical potentials. Then we have 

22A^/ + 13m^ ^' 22A^/ + 13m^ ^ ^ ^ 

On the other hand, in the broken phase (T < Tc), because of Q = and 
fit) = 0, we have, if the sphaleron process is in equilibrium, 

_ 8iV, + 4^ + 8 ^ _i6iV, + 9.n+18 

24Ar/ + 13m + 26 ^ ^ 24iV/ + 13m + 26 ^ ^ ^ ^ 

These equations imply that if the primordial B — L is absent, we have no 
baryon number and lepton number. Hence to have nonzero BAU after the 
EWPT, starting from the baryon-symmetric universe, 

(i) we must have B — L before the sphaleron process decouples, or 

(ii) B + L must be created at the first-order EWPT, and the sphaleron 
process must decouple immediately after that. 

Case (i) may be realized by mechanisms, stated in the previous section, such 
as GUTs, models with Majorana neutrinos, and the Affleck-Dine mechanism. 
The possibility (ii) is our main subject. 

It should be noted that to derive the above relations, we neglected the 
effects of mass in ( p.ll| ). In the standard model and its minimal super- 
symmetric extension, this mass effect can leave nonzero BAU starting with 
vanishing primordial B — L, if the model satisfies some constraints on the 



lepton-generation mixing or i?-parity violation [33 
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4 Electroweak phase transition 



Gauge symmetries broken by expectation values of scalars in nontrivial repre- 
sentations of the gauge group are known to be restored at high temperatures [0^ 
The nature of the phase transition depends on the parameters in the theory. 
Our aim is to clarify the properties of the EWPT, especially its dynamics. 
As is well known from the phase transition of classical systems such as spin 
systems and liquid-vapor phase transition, the dynamics of phase transitions 
are studied by use of the free energy as a function of temperature and or- 
der parameters. Fist we shall review analytic and numerical studies of the 
EWPT, which focus on equilibrium properties, and then some attempts to 
understand the dynamics of the EWPT. These are indispensable to determine 
whether the electroweak baryogenesis is possible and to estimate generated 
baryon number quantitatively. 

4.1 Static properties of the phase transition 

As noted in the previous section, the sphaleron transition in the symmetric 
phase and other elementary interactions are in equilibrium at T ~ Tc, be- 
cause of slow expansion rate of the universe. Hence static properties of the 
EWPT, such as the transition temperature, the order of the transition and 
the latent heat when the transition is first order, are described well by the 
equilibrium statistical mechanics. For this purpose, we calculate the effec- 
tive potential, which is the free energy density, as a function of the order 
parameters and temperature. In the case of a first-order phase transition, 
the effective potential has the form depicted in Fig. ^ in which 

y^c = lim y.(T) ^ (4.1) 

characterizes the order of the transition, with <f{T) being the order parame- 
ter. 

The first attempt to evaluate the effective potential was based on the 



loop expansion in finite-temperature field theory |35| , |36[| . The effective 
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Figure 2: Qualitative form of the effective potential at several temperatures. 
The potential at (/? = is normalized to be zero for all temperatures. 

potential in the minimal standard model (MSM), at the one- loop level, is 
given by 

Kff(<^;T) = \/trcc(^) + V^^\^-T), (4.2) 
where V|;ree(v') is the tree-level potential, 

Kree(^) = + (4-3) 

V^^^i^-^T) is the one-loop contribution 

V^'\^-T) = -WcAf logdet \w-/{k-^)\ (4.4) 

and the order parameter ip is introduced as the expectation value of the Higgs 
field as 

= ^ (° ) . (4.5) 

Here /ig and Aq are the bare parameters, which will be determined once 
one prescribes the renormalization. In( |4.4]) , the subscript A runs over all 
the species including the Faddeev-Popov ghosts, ca denotes the degrees of 
freedom of each particle and its statistics, that is, > for bosons and 
< for fermions, and is the inverse propagator, which is in general 
a matrix with Lorentz and internal symmetry indices and depends on ip 
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through its mass. The integration over the momentum should be understood 

as 

d^k 1,0 _ _ \ 2n7rT for bosons, 

(2^ with k = = I + i^^T for fermions. 

(4.6) 

For example, cw = 2 and iV-^'^{k;<^) = (-A;^ + m^(^))r/^'^ + (1 - ^1)^''^^ 
for W boson in the -Rg-gauge, where mw{'^) = \gVi c/ = —2 and iVj^{k; ip) = 
^ — mf{(f) for a Dirac fermion with mf{(f) = yf(f/^/2. The divergences in 
( |4.4| ) are absorbed in the bare parameters by the renormalization at zero 
temperature. For the Higgs boson, m\{Lp) = 3A(/?^ — /i^, so that V^^\ip]T) 
becomes complex for small v^[^]. Since the free energy is originally a real 
quantity, this pathology will be caused by the loop expansion, which is in- 
valid for small (p, and is expected to be cured if one takes the higher order 
contributions into account. At this point, we neglect the contributions from 
the Higgs and higher orders and discuss the structure of Kfr(v'; T). Now we 
consider only W and Z bosons and top quarks, since the contributions from 
the other fermions are negligible. Then the one-loop effective potential is 



V,f,{^;T) = Vo{^) + Vi^;T), 



(4.7) 



where 



rp4 

2^ 



[6lB{aw) + 3/B(az) - Glpiat)] ■ 



(4.8) 
(4.9) 



Here we follow the renormalization convention of Ref . |^ , and 

3 



B 

IbAo) 



{2m^ + mz- Am^] 



POO , 

J dxx^ log [1 =F 



where Vq = 246GeV is the minimum of Vo{ip) and = mA{'-p)/T. 



(4.10) 
(4.11) 
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For high temperatures T > mwjmz^mt, the integrals Ib,f can be ap- 
proximated by the high-temperature expansion! 



where 



V,s{^; T) ^ D{T' - T^)ip' - ET^p' + (4.12) 



D = A(2m2^ + m| + 2m2), (4.13) 



E = --3(2m3, + m|)~10-2, (4.14) 
1 



^0 

To = ^if^'-^Bv'o). (4.15) 

-4m^log^) (4.16) 

with loga^ = 2 log An — 2je, loga^? = 2 log n — 2'~^e and 7^; = 0.5772- ■ • 
being the Euler constant. The effective potentials of other models which ex- 
hibit first-order phase transition will have a form similar to that of ( [4.12| ) with 
-E > 0, while each coefficient will be modified appropriately. By use of this 
approximate form, we see that aX T = Tc ■, where the transition temperature 
Tc is determined by 

the effective potential has two degenerate minima at (f = and 

VC = ^. (4,18) 

The presence of the y^^-term with positive E in the effective potential is 
essential for the EWPT to be first order. It originates from the zero-frequency 
mode of the bosons in the finite-temperature Feynman integral in ( [4.9|) . This 



can be seen as follows p7|. The zero-frequency contribution to dVes/dcp has 
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the form of 



dip dip J k'^ + 171"^ (p) An dip 



VVn%), (4.19) 



where the divergence which is absorbed by renormahzation is dropped. Since 
this is proportional to p"^, V^s contains a p^-term upon integration. Note that 
within this approximation, the EWPT appears to be first order, since Ves 
always has a y^^-term whose coefficient has the correct sign. 

Once the expectation value of the Higgs at Tc is known, one can evaluate 
the sphaleron rate at Tq- Requiring that the sphaleron process decouples 
in the broken phase just after the EWPT imposes an upper bound on the 
Higgs mass[^]. This amounts to the condition that the sphaleron rate in the 
broken phase ( |3.7D multiplied by particle density ~ is smaller than the 
Hubble parameter (|3.9|) , which yields 

g > 1- (4.20) 

This is converted to an upper bound on A. Since the Higgs mass is given by 
iriH = V^Xvq at the tree level, rriH is bounded as 

mn < 46GeV. (4.21) 

This is obviously inconsistent with the present lower bound mn > 58.4GeVp|]. 
Although this upper bound may be a crude value, the decoupling of sphaleron 
process in general will require the Higgs be lighter than some bound. This 
upper bound will be relaxed if one extends the model to include extra bosons, 
which effectively enhance the p^-teim. 

Among these extended models are those with two Higgs doublets, includ- 
ing the minimal supersymmetric standard model (MSSM). We studied the 
massless two-Higgs-doublet model, in which the gauge symmetry is broken by 
the Coleman- Weinberg mechanism, without recourse to the high-temperature 
expansion. Our result shows that for the entire range of the parameters 
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we considered, the EWPT is first order, the sphaleron decouphng condi- 
tion is satisfied and the hghtest Higgs scalar is as heavy as 63GeV[^. The 
more general two-Higgs-doublet model was studied by randomly scanning 
the vast parameter space |HD|. For the parameters consistent with present 
experiments and allowing perturbation, the lightest Higgs scalar has mass 
45 < TUh < 190GeV. The upper bound on the Higgs mass in the MSSM 
was evaluated for various values of the soft supersymmetry-breaking param- 
eters and was found to be near the experimental lower bound|^. We note 
that all these results are based on the one-loop effective potential and the 
latter two works used the high-temperature expansion. In general, for these 
extended models, we have more chances to satisfy the bound on the Higgs 
mass consistent with experiments and to obtain first-order EWPT. 

Here we shall comment on higher-order studies of the effective poten- 
tial. The most dominant contributions are the daisy type diagrams, and 
summing up these gives a well-defined effective potential ||3^. Since these di- 
agrams correspond to infinite-order mass corrections, some authors evaluate 
the 'improved' effective potential simply by replacing of the particle in 
the loop with + Il{(p) in the one-loop effective potential, where is 
the one-loop polarization. Although this prescription may somehow cure the 
pathology of complex effective potential for smaller (p, it does not yield correct 
answer [p7| . The correct method gives a smaller value of E in ( [4.12| ), which 
means a weaker first-order phase transition. Some two-loop calculations sug- 
gest the opposite situation. The two-loop corrections to the ring-improved 
one-loop effective potential were found to raise the lower bound derived from 
the sphaleron decoupling condition, while it still rules out the MSM||4^.^ 
The two-loop effective potential was shown to increase the strength of the 
first-order EWPTp4|. These two-loop corrections are so large that we still 
need more reliable evaluation for a definitive conclusion. 



*This calculation was carried out in the Landau gauge. As for the gauge-dependence 



of the effective potential, see Ref. [43 
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Besides these analytic investigations, the static properties of the EWPT 
have been studied numerically by the lattice simulations [^]. Since at present 
the electroweak theory cannot be put on the lattice completely, its subsys- 
tem, SU{2) gauge-Higgs system with one Higgs doublet, has been taken, 
and two methods have been applied to the simulation. At high tempera- 
tures, fermions attain effective masses of 0{T) because of their antiperiod- 
icity along the imaginary time direction, even if they are massless at zero 
temperature. This fact might justify this approximation neglecting all the 
fermions in the model at high temperatures. One of the simulations is the 
standard one on the periodic euclidean four-dimensional lattice. The other 
is that of the effective three-dimensional theory, which is obtained by taking 
the high-temperature limit. In this limit, the euclidean time period, 1/T, 
is smaller and the system is reduced to that of the fields, which are zero 
modes of the Matsubara frequency in the original system. That is, the ef- 
fective three-dimensional system is composed of an SU{2) gauge field with 
one Higgs doublet and one Higgs triplet, which is the remnant of the time- 
component of the original gauge field. Correspondingly, the parameters come 
to have temperature-dependence, which are determined by evaluating finite 
temperature Green's functions in both the original and the reduced theory 



and by matching them[^]. Because of Elitzur's theorem, no local symmetry 
is broken on the lattice, so that (0) = 0[^. The measured quantities are 
expectation values of the gauge invariant operators such as {4>^4>), its jump, 
transition temperature Tc, latent heat and surface tension. The latent heat 
is obtained by the derivative with respect to temperature of the difference of 
the free energy between the two phases, just as the continuum theory. There 
are three methods to measure the surface tension between the two phases. 
As for the comparison of the results of these simulations with the continuum 



analysis, see Ref. [|^. According to Ref. [Q, the numerical results coin- 
cide with those of the continuum two-loop perturbation theory, up to the 
Higgs masses of about 70GeV. For example, the three-dimensional reduced 
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theory gives, with both perturbative and numerical methods, Tc — 93GeV 
and i^c/Tc ^ 1.8 for m// = 35GeV, Tc ~ 140GeV and i^c/Tc ^0.7 for 
mn = 60GeV, Tc ~ 155GeV and (fic/Tc ^ 0.6 for mn = 70GeV. It also 
yields the result that if we denote the surface tension by a, <j/T^ ~ 0.84 
for niH = 18GeV, a/T^ ~ 0.065 for uih = 35GeV and a/T^ ^ 0.008 for 
mn = 49GeV by the four-dimensional simulation. From these results, the 
EWPT of the MSM is first order for m/^ < 70GeV. The strength of the 
transition rapidly decreases as mn increases. For ttlh > 60GeV, ( [4.20|) is not 
satisfied, but the upper bound on ttih might be weakened since temperature 
in the broken phase is lower than Tc because of supercooling. 

These quantities compared between analytic and numerical methods are 
those derived from the effective potential but not itself. One can, in principle, 
calculate the effective potential by numerical methods, but its form is always 
convex p6[| so that it does not give quantities such as the height and width 



of the barrier of the effective potential dividing the two phases. Our aim is 
to obtain a 'classical' effective potential which offers information about the 
first-order phase transition. Such an effective potential should reproduce the 
results of the lattice simulations and have a shape like Fig. ^ 

As we shall see in the next section, to have an efficient CP violation, 
extension of the Higgs sector would be needed. Further, such extension would 
open the chance to fulfill the condition ( [4. 201) within the present Higgs mass 
bound. Hence both analytic and numerical studies of such models would be 
desired. 



4.2 Dynamics of the phase transition 

To have nonzero BAU, we need a nonequilibrium state, which is realized by 
expanding bubbles created at the first-order EWPT. Whether the EWPT 
is a first-order transition accompanying nucleation and consecutive growth 
of bubbles of the broken phase in the symmetric phase will be determined 
by the form of the effective potential. If the EWPT proceeds with bubble 
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nucleation and growth, the width and velocity of the bubble wall are essential 
to evaluate the baryon number generated. The nucleation rate of the bubbles 
is also crucial, since if their nucleation dominates over their growth, the total 
region which experiences nonequilibrium conditions will decrease. 

A crude picture of how the EWPT proceeds can be gained as follows. 
Suppose that the effective potential at temperatures around Tc is known. 
The nucleation rate per unit time and volume is given byp6 



/(T) = /oe-^^(^)/^, (4.22) 

where the prefactor Jq is determined by temperature, viscosity and correlation 
length in the symmetric phase and AF is the change in the free energy caused 
by formation of a bubble. As long as the thin-wall approximation is valid, 
AF is approximately given by 

A-jr 

AF{T) = —r^[ps{T) - Ph{T)] + Anr'a, (4.23) 

where ps{b) is the pressure in the symmetric (broken) phase and is given by 

Ps{T) = -Kff(0;T), p,(T) = -V;ff(^(T);T). (4.24) 

In general, Ps(T) < PhiT) because of supercooling, a is the surface energy 
density given by a ~ / dz{dip/ dz)"^, with z being the coordinate perpendicular 
to the bubble wall. For smaller radius r, the second term on the right-hand 
side of ( [4. 231) dominates so that the surface tension shrinks the bubbles. For 



larger r, the first term dominates and the bubble expands to lower the free 
energy. The critical bubble has a boundary radius between these bubbles, 
whose radius is 

r.(T) = (4.25) 

With the free energy Ves{<^',T) known, one can calculate various thermo- 
dynamical quantities. For example, the entropy density is given by s = 
—dVcs/dT and the energy density is p = Ves — Ts. The latent heat is the 
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difference of the energy density between the two phases. How the phase 
transition proceeds may be characterized by the fraction of space which has 
been converted to the broken phase, f(t). It obeys the integral equation 

fit) = f dt' I{T{t'))[l - fit')]V{t',t), (4.26) 

Jtc 

where V{t',t) is the volume of a bubble at t which was nucleated at t' < t. 
Here, the time t and temperature T are related by the scale factor of the 
universe R at t, since p oc for the radiation-dominated universe and 
p = {tt"^ / 30) g^T*. (See Appendix A.) V(t',t) is taken as 



V{t',t) g 



n{T{t'))+ / dfv{T{t")) 



n 3 



(4.27) 



where v{T) is the wall velocity. This simple approach does not take into 
account interactions and fluctuations of the bubbles. The wall velocity 
and thickness are estimated by solving dynamical equation for (f{x) with 
friction [49, 5C]. In Ref.||5l[], ( [4.26 ) is numerically analyzed by use of the wall 
velocity given in Ref. |49| with the one- loop improved effective potential of 
the MSM for mn = 60GeV and rrit = 120GeV. The results show that if we 
measure the time after the temperature reached Tc, (1) at 6.5 x 10~^'^sec, 
bubbles began to nucleate,^] (2) at 6.87 x lO^^^sec, the nucleation was turned 
off and then only about 10% of the universe had been converted to the broken 
phase, then (3) the remaining 90% of the universe is converted by the bubble 
growth with almost constant velocity of about 0.8, and (4) the transition is 
completed at 7.05 x 10~^^sec. Further if we denote the temperature at which 
the bubble nucleation begins by T/v, (Tc — Tn)/Tc — 2.5 x 10~^, which means 
the EWPT shows small supercooling. 

All these results are based on the effective potential of the MSM at the 
one- loop level or the improved version of it. As mentioned, the above esti- 
mation neglects the interactions of the bubbles and various fluctuations. For 

^At this temperature, the characteristic time scale of the electroweak processes is 
O(10-26)sec. 
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weak first-order transitions, thermal fluctuations might affect the dynamics 
of the transition. Then the barrier dividing the symmetric and broken phases 
of Ves is so low that the two phases may be mixed by thermal fluctuations. 
In such the EWPT might proceed not with nucleation of the criti- 

cal bubbles but with that of the subcritical bubbles[^] so that the simple 
analysis above may no longer be valid. 

The numerical simulations as well as the two-loop effective potential sug- 
gest a stronger flrst-order transition in the MSM than the one- loop result. 
Extensions of the standard model with more Higgs scalars may also lead to a 
stronger transition, since the upper bound on the Higgs mass is raised some- 
how and since this also implies that (pc is larger for a given Higgs mass. In 
fact, our work on the two-Higgs-doublet model supports this fact and thinner 
bubble walls . 



5 Mechanism of the electroweak baryogene- 
sis 

If the sphaleron process is in equilibrium in the symmetric phase and de- 
couples in the broken phase at the EWPT, it is out of equilibrium around 
the moving bubble wall dividing the two phases. In this section, we review 
how the baryon number is generated by use of this nonequilibrium situation. 
Around the expanding bubble wall, the expectation value of the Higgs fleld 
is schematically depicted in Fig. |^, where ip represents the absolute value of 
the Higgs scalar, vq is that in the broken phase and Vco is the value above 
which the sphaleron process decouples. There are two mechanisms of the 
electroweak baryogenesis; classical or adiabatic mechanism, which is called 
'spontaneous baryogenesis' and quantum mechanical one, which is also called 
'charge transport scenario'. Although these two mechanisms may be respon- 
sible for baryogenesis at the same time, the former is effective for a thick 
bubble wall while the latter is effective for a thin wall, as we shall see be- 
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Figure 3: Profile of tlie expanding bubble wall. In the region right to the 
vertical dashed line, the sphaleron process is effective, is the wall velocity. 

low. Further both require extensions of the Higgs sector in the MSM, except 
for the scenario proposed by Farrar and Shaposhnikov[^]. The spontaneous 
scenario was proposed earlier than the charge transport mechanism. It, how- 
ever, was pointed out that it would generate too small a number of baryons 
in its original form. After the charge transport scenario was proposed, the 
classical mechanism was revived by introducing a nonlocal effect, which is an 
advantage of the charge transport scenario. 

5.1 Charge transport scenario 

Fermions interact with the bubble wall through the Yukawa couplings, which 
in general contain CP violation. If the CP violation is not constant over 
spacetime, it cannot be rotated out by a biunitary transformation so that 
there remains a physical effect. This CP violation discriminates between 
the interaction of the fermions and that of the antifermions with the bubble 
wall. That is, the reflection rate of the fermions does not coincide with that 
of antifermions so that the net current of some quantum number flows into 
the symmetric phase region, where the sphaleron process is effective. There 
the previous equilibrium state will be forced to shift to new state, which will 
contain baryon number excess. One can see how this occurs by solving the 
kinetic equations, which will contain various interactions of different time 
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Figure 4: The free energy as a function of the baryon number is biased and 
a nonequihbrium stationary state is reahzed. 



scales. Such equations are in general difficult to solve. When the bubble 
wall moves with constant velocity, the flux will bias the free energy along 
the baryon number to realize a stationary nonequihbrium state, as depicted 
in Fig. 1^. In this case, the constant flow of the flux will induce a nonzero 
chemical potential for the baryon number. Then the steady state produces 
the baryon number according to the equation 

UB - (5.1) 

where Fsph is the sphaleron rate in the symmetric phase (|3.8| ) . For the deriva- 
tion of this equation, see Appendix B. 

Cohen, Kaplan and Nelson related the flux flowing into the symmet- 
ric phase to the generated baryon number as follows[Q. Suppose that a 
left(right)-handed fermion of species i with charge Q^l(^^^ is incident from the 
symmetric phase region. The expectation value of the injected charge in the 
symmetric phase, which has been brought by the reflection and transmission 
of the fermions, is given by 

AQI = [{Q'r-QI)R1-,r + {-Q1 + Q'r)RI-,r 

H-Q\){Tl^L + ti^r) - i-QR)in^L + n^R)]fL 
+[iQl - Qr)Rr^l + {-Qr + Q'l)RUr 
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H-Q'r)Tr^l + TI,^r) - {-Qim^L + n^R)]fm (5.2) 

where Rp^^i {Rfi^l) is the reflection coefficient for the right-handed fermion 
(antifermion) incident from the symmetric phase region and reflected as 
the left-handed one, T's are the transmission coefficients and is the 

Ieft(right)-handed fermion flux density in the symmetric phase, which will 
be given below. Here we have used f[ = In a similar manner, the expec- 
tation value of the change of the charge brought by the transmission of the 
fermions incident from the broken phase region is 

^Q'^ = QUTLLfL + T'n^LfR^+QRiTLRfL + T'n^RfR^ (5.3) 
H-QDin-^Lfl + TLlIr^ + i-QR)in-.RfL + T^RfRi)- 

By use of the unitarity 

Rl^n + T[^L + Tl^R = h etc. (5.4) 



and the reciprocity ||55|| 



n^L + n-.R = TLL + T'n^L, etc. (5.5) 
the total expectation value of the change is 

AQt + /\Q\ = {Ql - Q\){n - f')^R, (5.6) 

where 

AR ^ R%_^^ - R%^^, (5.7) 

and we have put //^^^ = //^^^ = Z/*'^''. The total charge injected into the 
symmetric phase can be obtained by integrating this in the rest frame of the 
bubble wall and going back to the rest frame of the medium, 

= %^ r r ^p^p^ [f'(PL,PT) - f!:{-PL,PT)] Ai?(^, ^). 

(5.8) 
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Here 

f!{PL,PT) 



PL 1 



m-PL,PT) = f ^ , (5.9) 



are the fermion flux densities in the symmetric and broken phase, respec- 
tively, mo is the fermion mass in the broken phase, is the wall velocity, 
7 = 1/ ^1 — v^, E = \Ip'l+ Pti Pt being the transverse momentum with 
respect to the bubble wall. If we calculate AR with zero-temperature Dirac 
equation in the presence of the CP-violating potential, it is a function of 
mo/ct and Pi/a, where is the wall width. 

Ai? contains the information of CP violation for the model with which we 
are concerned. In the MSM, CP violation enters only through the Kobayashi- 
Maskawa (KM) matrix in the charged-current interactions. Farrar and Sha- 
poshnikov proposed that even in the MSM, the interactions of the quarks with 
the plasma of the gauge bosons and Higgs bosons, which modify the disper- 
sion relations of the quasiparticles, may yield significant baryon number|]53[ . 
The effect of the KM matrix comes into the dispersion relations of 0{aw) and 
there are claims that QCD damping effects, whose scale is shorter than that 



of other interactions, may bring quantum decoherence|[56||. Hence to obtain 
enough BAU only via electroweak baryogenesis, we would have to employ 
an extra source of CP violation, such as that in the Higgs sector in the 
two-Higgs-double model. The CP violation in the Higgs sector affects the 
fermion propagation around the bubble wall at the tree level, so it was the 



first model considered in the charge transport scenario|^J]. As for the quan- 
tum number injected in the symmetric phase, we must take that conserved 
there and satisfying Ql — Qr 0. One such charge is the weak hypercharge 
Y . Now let us see how the injected hypercharge biases baryon number. 

As we stated above, the change of state caused by the injection of hyper- 
charge would be described by the kinetic equations. Here we assume that the 
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bubble wall is macroscopic and moves with almost constant velocity, and that 
the state shifts to a new steady state. Such a situation would be described 
by the chemical potentials for all the particles in equilibrium. This simple 
picture would be valid in the region far from the wall and in the case that 
the elementary processes are fast enough to realize a new stationary state. 
For simplicity, we consider only the third generation and the case that all 
the particles had been in equilibrium with ub — riL — through the gauge 
and Yukawa interactions before the injection of the hypercharge. Now we 
introduce the chemical potentials 

l^tL^l^bL: l^tR: l^bR: I^tl : l^>^r ^ I^tr, I^'W for W', /Xq for 0°, /X_ for (f)_. 

When the charged-current interaction is in chemical equilibrium, 

and when the Yukawa interaction is in chemical equilibrium, 

1^0 = -l^ti^ + l^tR = -l^br. + l^bR = -l^rr. + I^tr- (5.11) 

No further new equations arc derived if the other elementary processes are in 
equilibrium. We also introduce the chemical potentials for quantum numbers 
which are conserved or almost conserved because of the slow processes: 

IJ'B-L, fJ'Y, fJ'Is', fJ'B- 

Here we assume that the sphaleron process is out of equilibrium. Otherwise, 
/iB vanishes so that the baryon number is not generated. The chemical 
potentials of the particles are expressed by these as 

11 1 , ,1 

l^tUbL) = :^I^B + -I^B-L + -I^Y + 

1 1 2 

1 1 1 

fJ'bR — i:^B + -^IJ'B-L - -^H'Y, 
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(5.12) 



Because of ( p.ll| ), the baryon and lepton number densities are given by 

ub = 3 ■ —(/ii^ + /ii^ + /ifo^ + /ife^) = —(2;Ub + 2;UB_i + /iy), (5. 13) 
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If?7,B = ''^L = before the injection of the hypercharge flux, 



^J^B-L 



_ 1 

— -fj'Y, 

o 



(5.14) 



(5.15) 



hold because of ( |5.13| ) and ( p.l4| ). Then the hypercharge density, which is 
expressed by the chemical potentials of the particles, can be given, with the 
help of (|5l2|) and {^J^, by 



Y 
~2 



( \1 2 1 , 

-^(m+-)/.y. 



^(/^i.. +/^ri) -/irflj 



(5.16) 



where m is the number of the Higgs doublets. From ( ^.15| ) and ( p. 16 ), we 
have 

Y 

= 2(m + 5/3)T2- ^^-^^^ 
Thus the injected hypercharge biases the free energy along the direction of 
the baryon number. By use of ( p.l| ), we have the generated baryon number 
density 



riB 



sph 



T 



dt jiB 



sph 



2(m + 5/3)T3 



dtY. 



(5.18) 
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If we denote the hypercharge density at distance z from the bubble wall by 
Py{z)., the integral of Y is estimated as 

/r-zjv^ \ poo 

dtY = dt py{z - v^t) = — dzpviz). (5.19) 

J-oo Vw JO 

Since the hypercharge supplied from the wall would reach a finite distance 
when it is caught up with the wall, the last integral above will be approxi- 
mated as 

1 r°° FvT 

— dzpriz)^^, (5.20) 

where r is the transport time within which the scattered fermions are cap- 
tured by the wall. Hence the generated baryon asymmetry is 

^^Arif ..a^^.4,.rT, (5.21) 

where A/" is a model-dependent constant of 0(1). r may be approximated by 
the mean free time or diffusion time D/vw with the diffusion constant D of 
the charge carrier, where D ~ 1/T for quarks and D ~ (10^ ~ 10'^)/T for 
leptonsfS^. It was shown that if the forward scattering is enhanced, even for 
the top quark, tT ~ 10 ~ 10'^ depending on f^, where the maximum value 
is realized at Vyj ~ l/-\/3|^. For this optimal case, 

^-10-»-45- (8.22) 

Then the hypercharge flux Fy/{vwT^) ~ O(10~^) would be sufficient to ex- 
plain the present BAU. For leptons, tT is much larger than that of quarks, 
since its propagation is not disturbed by strong interactions. As we shall see 
below, the hypercharge flux depends not only on the mass of the carrier but 
also the wall width and velocity. 

Now we evaluate the chiral charge flux in ( p.8|) . Here we assume that 
the simple zero-temperature Dirac equation describes the propagation of 
fermions well. This is justified when the mean free path of the fermion is 
larger than the bubble wall width. Since the dynamics of the EWPT in the 
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extended models is not well established, we treat the wall width and velocity 
as parameters and calculate the flux as a function of them. The problem 
is to solve the Dirac equation in the background of the bubble wall. The 
bubble wall is composed of classical gauge and Higgs fields. As we shall see 
in the next section, its profile will be composed of the Higgs scalar only and 
the gauge fields are pure-gauge type so that we fix the gauge such that the 
gauge fields vanish. Then the Dirac equation to be solved is, for one flavor. 



i0ip{x) — m{x)PR%lj{x) — m*{x)PLip{x) = 0, 



(5.23) 



where Pr = ^-^^Pl = and m{x) = — /(0(x)) is a complex- valued func- 
tion of spacetime, with / being the Yukawa coupling. When the radius of 
the bubble is macroscopic and the bubble is static or moving with a con- 
stant velocity, we can regard m(x) as a static function of only one spatial 
coordinate: 

m{x) = m(t, x) = m{z). 

Putting 



^(x) = ei'^i-Et+PT-XT) ^ = +1 or - 1) 



( |5.23| ) is reduced to 

(t(7°£' - 7tPt) + il^dz - fnn{z) - 2757717(2;)] iPe{Pt, z) = 0, 
where 



(5.24) 



(5.25) 



Pt 

ItPt 
miz) 



(p\p2^ 



Xt = {x ,x), pt = \Pt\ 



tV^ + 7V, 



If we denote E = E* cosh 77 and px = E* sinh7/ with E* = y E'^ — pj^, px in 
( |5.25| ) can be eliminated by the Lorentz transformation ip ip' = e~^'^°'^'^'ip: 



dzipEiz) 



77 



(tE*7° + mFi{z) + i'y^miiz) ipEiz). 



(5.26) 
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Figure 5: AR for the potential (5-27). We take mo = 1. 

This equation was first solved numerically in a finite region containing the 
bubble wall |^4|| , with the following potential 



, , 1 + tanh(az) / 1— tanh(az)\ 
m{z) = mo ^^-^ exp ( -ztt 1 • (^-27) 

The nontrivial 2;-dependence of the phase of m{z) yields CP violation in the 
Dirac equation. If it were independent of z, the phase could be eliminated by 
redefinition of the wave functions, so that it would yield no physical effect. 
We analyzed ( |5.26| ) by the perturbative method]^, which can be applied 



when the imaginary part of the mass function is much less than mo, as well 
as by numerical methods [Q, which solves ( |5.26| ) in the infinite region and 



gives the results with any precision. In Fig. ^, we show data for AR with the 
potential (|5.27|) , which were obtained numerically . For higher incident 



energy, |Ai?| decreases exponentially. This is because for pi > mo, the 
wave function is obtained approximately by the WKB method, which infers 
that Rr^l and Rr-*l are exponentially small. \AR\ can be 0(1) when the 
wall width is comparable to the Compton wave length of the carrier, i.e., 
mo /a ~ 0(1). Note that larger Yukawa coupling does not always yield 
larger asymmetry in the reflection coefficients. The chiral charge flux, which 
is essential to determine the generated BAU, is evaluated by inserting AR 



34 



logio(m/r) 
2 




logm(a/T) 

Figure 6: Contour plot of the chiral charge flux, normalized as 
^oSio [~F'qI {v^T'^iQi — Qr))] for the potential (5-27). Here we take v^, = 0.1 
and T = lOOGeV. 



into ( |5.8| ). We normalize it in a dimensionless form as 

Fo 



(5.28) 



uT^iQL-QRV 

which enters in the baryon asymmetry of (|5.21|) or ( |5.22| ). The numerical 
results for the potential (|5.27| ) are shown in Fig. Fig. ^ and Fig. | for the 
wall velocities = 0.1, = l/VS ~ 0.58 and = 0.98, respectively, as 
functions of the carrier mass and the wall width. Here we take T = lOOGeV. 

These figures suggest that the flux rapidly decreases for a heavier carrier 
and thicker wall and the region in which —Fq/{vwT^{Ql — Qr)) > 10~^ 
becomes broader for larger wall velocity. The maxima in these figures are 
realized at a ~ T, for which the wave length of the carrier is comparable 
to the wall thickness. For a thick bubble wall a ~ 40/T, fermions of mass 
O.IT < tuq < T can yield a flux larger than 10~^, while for a thin wall case 
a ^ 1/T, fermions of mass 0.03T < mo < lOT can yield sufficient flux. 

To calculate the flux, we assumed the form of the complex mass as ( ^.27 ). 
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Figure 7: Contour plot of the chiral charge flux, normalized 
^^Sio [~Fq/ iy^T'^iQi — Qr))] for the potential (5-27). Here we take 
0.58 and T = lOOGeV. 
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Figure 8: Contour plot of the chiral charge flux, normalized 
logio \~PqI{^w^'^{Ql — <5fl))] the potential (5-27). Here we take v, 
0.98 and T = lOOGeV. 
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It yields no CP violation in the broken phase so as to be free from any 
experimental limitations. However, the profile of the bubble wall should be 
determined dynamically. This is the subject of the next section, and it will 
also affect the other scenario of the electroweak baryogenesis. 

5.2 Spontaneous baryogenesis 

This scenario, in contrast to that in the previous subsection, is classical and 
effective for thick bubble walls. Although this scenario has several variants, 
the first proposed may be that by Shaposhnikov[^], That case depends on 
the existence of the condensation of the Chern-Simons number, which adds 
a CP- violating term in the effective action as 

AS'eff = ^^Ncs. (5.29) 

where Ncs is defined by ( |3.4| ). This scenario fails in the MSM since it can 
induce at most /i = 10^^^ while /i = 10"'' is needed to explain the BAU 



It is later elaborated by noting that in the two-Higgs-doublet extension of 
the MSM, the effective action takes the form as 

A£eff = -^e{x)F,,{x)F'^^{x), (5.30) 

where 6{x) is the relative phase of the two Higgs doublets [^, This term is 
CP- invariant since d{x) is CP odd, so that it is free from any constraint 
on CP violation. A nontrivial evolution of 6{x) during the EWPT would 
induce asymmetry like /i in ( p.29| ), upon integrated by parts. For high tem- 



peratures, the factor in front of the Chern-Simons number is accompanied 
by {rrit/TY coming from the top quark loop[^. Cohen et al. proposed a 
different approach, which utilizes the bias for the hypercharge in place of the 



Chern-Simons number, in the two-Higgs-doublet extension of the MSM|]6T 
Their observation was that if the neutral component of the two Higgs scalar 
takes the form 

0°(x)=p,(x)e^^^ (j = l,2) (5.31) 
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where only 0i is supposed to couple to the fermions, one can eliminate 6i from 
the Yukawa coupling by an anomaly-free {U{1)y) rotation, which changes the 
fermion kinetic terms by 



1 _ 2 _ 1 - 

-qL{x)YqL{x) + ■::Un{x)YuR{x) - -dii{x)YdR{x) 

6 6 

—hix^'^lLix) - enix)YeRix)^ , (5.32) 



where summation over generation is implicit. During the EWPT, 6i{x) will 
behave nontrivially so that (^i) ^ 0. Then (^^i) will play the role of 'charge 
potential'. If this biases the hypercharge in the symmetric phase region, the 
same procedure to generate the baryon number as discussed in the charge 
transport scenario works here. They estimated the generated baryon asym- 
metry as 

^ ~ 10-^ A9, (5.33) 
s 

where A6 is the change of 6i during the EWPT. None of these mechanisms 
rely on any quantum mechanical scattering, so that their effects are not 
weakened for thick bubble walls, in contrast to the charge transport scenario. 

It was pointed out that these scenarioes cannot generate sufficient baryon 
asymmetry, because of an additional suppression factor in the charge potential |62 
Since in the symmetric phase the CP violating angle 6i does not affect the 
dynamics because of pj = 0, the charge potential will be zero there. The 



hypercharge current in ( |5.32|) is the fermionic one, and taking correctly the 
contribution from the scalar parts, the generated term by the rotation is 
dfiOijy with the total hypercharge current jy, which is conserved in the sym- 
metric phase. On the other hand, its nonconservation in the broken phase 
leads to, at high temperatures, 

2 

d,d,-jl^^'^F,,F^\ (5.34) 

However, the sphaleron process is enhanced in the symmetric phase and in a 
restricted region in the broken phase where the expectation value of the Higgs 
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is so small that the sphaleron rate is large. In Fig. such a region is that 
where v < Vco- Hence the estimate of the generated baryon number suffers 
from a suppression factor {vco/TY, which could be as small as 0(10^^). 

Later Cohen et al. found that the effects of diffusion can enhance the 
baryon asymmetry |]ti3||. The diffusion may carry the charge asymmetry into 
deep symmetric phase region. They solved the diffusion equations for various 
charges, assuming the profile 



1 — tanh(az) 
\(p\z)) = V exp 



. . ^ 1 — tanhfa^;) 
-iA9 ^ — - 



(5.35) 



with A6 = 7r/2, and found that the generated baryon asymmetry becomes 
enhanced by a factor of 0(1/ a^) ~ 10^ over the previous estimates, and that 
the baryon asymmetry is insensitive to the details of the bubble dynamics, 
such as the wall width 1/a and velocity Vw Once diffusion is taken into 
account, generation of the baryon number is similar to that in the charge 
transport mechanism. In this sense, these scenarios are also called 'nonlocal 
baryogenesis'. For recent and more elaborated studies on these subjects. 



which include the effects of diffusion, see Ref.[§4 . 

As we saw in this section, to obtain enough BAU, some extensions of 
the MSM would be needed, in both mechanisms. How these mechanisms 
contribute to generate the BAU would be determined by the dynamics of the 
EWPT. We present some results based on the charge transport mechanism. 
For thicker bubble walls (smaller a/T), the other mechanism will come into 
effect, so that the flux for smaller a/T may not damp so rapidly in practice. 
In the illustrative calculations above, the profiles of the CP violating bubble 
wall, ( |5.27D and ( ^.35| ), were assumed without any reasoning. They should 



be determined by the dynamics. We need to relate CP violation around 
the bubble wall to the parameters in the models of electroweak theory, to 
estimate quantitatively the generated baryon asymmetry. 
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6 CP- violating bubble wall 

6.1 The model and the equations of motion 

C and CP violation are among the requirements to generate baryon asym- 
metry starting from a baryon-symmetric universe. C is violated in any chiral 
gauge theory such as the standard model, while CP is violated only through 
the KM matrix in the MSM. As we saw in the previous section, some ex- 
tension of the Higgs sector of the MSM would be needed to have a sufficient 
baryon asymmetry only by the electroweak baryogenesis. In practice, the 
amount that the baryon number is generated depends on the CP violation 
of the model under consideration. As we noted, one of the most attractive 
features of the electroweak baryogenesis is that it relies exclusively on physics 
which can be verified by experiments in the near future. Our aim is to de- 
velop a formalism by which the generated baryon number can be evaluated 
quantitatively once one selects the model lagrangian for the electroweak the- 
ory. Here we concentrate on the two-Higgs-doublet models, which contain 
the minimal supersymmetric standard model (MSSM) as a special case and is 
the simplest model to incorporate the Higgs-sector CP violation. Assuming 
the effective potential which gives the first-order EWPT, we determine the 
functional form of the CP violating phase in the Higgs fields. For a review 
of this model and MSSM, see Ref. . 

The most general renormalizable Higgs potential which is invariant under 
the SU{2)l X U{1)y gauge transformation is 

Vq = m^$|$i + ml^l^2 + {ml^\^2 + h.c.) 



where the charge assignment of the scalars are (2,1) for SU{2)l x U{1)y- 
Here the hermiticity of Vq requires ml, ml, Ai, A2, A3, A4 G R while in general 
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"^3' -^5) -^6) -^7 £ C, three of their phases are independent and yield the exphcit 
CP violation. The Yukawa interactions are generally 

+rL(/f^$i + #$2)efi + h.c., (6.2) 

where I'j = ia2^* is the charge-conjugated Higgs fields and fi^^f'^^ are the 
Yukawa couplings which are matrices with the generation indices. If we 
impose the discrete symmetry, 

ur ^ -ur, (Ir ^ dR, cr ^ cr (6.3) 

or 

Ur ^ Ur, dR dR, Cr y-^ Cr (6.4) 

which forbids the tree-level Higgs-mediated fiavor-changing neutral current 
(FCNC) interactions ||66|. the Yukawa interactions in ( |6.2D are restricted to 



the form of 

Cy = QLl^'^^^idR + qL&^^2UR + iLf'^^^ieR + h.c. (6.5) 

or 

Cy = qLf^^'^^idR + qL&'^^iUR + hf^'^^ieR + h.c, (6.6) 

and at the same time, in the Higgs potential, it is required Ag = A7 = m\ = 0. 
One can introduce m| 7^ 0, which breaks the discrete symmetry softly, since 
it does not induce divergent counterterms of the Aej-type so that there is 
still no tree-level FCNC interactions. 

In the MSSM, the tree-level Higgs potential is given by ( |6.1|) with 

Ai = A2 = \{g' + g"), A3 = \ig'- g"), X, = \g\ 

A5 = As = A7 = 0, (6.7) 
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and the Yukawa coupling is of the type (|6.5|), where is identified with $i 
and $n with $2- The m|-term in the MSSM comes from the soft supersym- 
metry breaking term and is generally complex. 

When ml, A5, Ag, A7 G R, which is the case for the MSSM, CP is not 
violated by the potential. However, CP can be spontaneously broken if 



A. < 0, 



2ml 



< 1, 



2X5V1V2 

where we parameterize the vacuum expectation value of the Higgs as 

\ 



(6.8) 



V2 



Vie 



9i 



72- 



(6.9) 



This is because, upon inserting 
becomes 

cos6' + 



Vo = -V^V2 



9|) into (|6.1|), the 6'-dependent part of Vq 



Ayf 2 — 2ml 



+ 



(6.10) 



2X5V1V2 

which takes the minimum at 9 ^ mr {n E Z) if the condition (|6.8| ) is 
satisfied.^ Although the MSSM does not admit the tree-level spontaneous 
CP-violation, the radiative corrections can induce the four-point couplings in 
such a way that ( |6.8|) is fulfilled ||6^. When some discrete symmetry is broken 
radiatively, one inevitably gets a pseudo-Goldstone boson, which is usually 
light because it acquires its mass only through radiative corrections! 70 1. In 
the case of spontaneous CP violation in the MSSM, the light scalar mass is 
about 6GeV at most, which is phenomenologically forbidden . 



We are interested in CP violation at finite temperature, especially near 
the bubble wall created at the EWPT. The effective potential of the model at 
finite temperature, which depends on the expectation values of all the Higgs 
fields, would determine to what extent CP is violated at T ^ Tc in the broken 
phase region. In fact, the one-loop calculation in the MSSM shows that 
CP could be spontaneously violated at high temperature while unbroken at 



^As for a solution to the domain wall problem associated with spontaneous CP viola- 
tion, see Ref. [pSl. 
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zero temperature, avoiding the light pseudo-Goldstone boson|7T|, |75|, |7^. We 
expect that if we know the global structure of the effective potential, we could 
derive from it the profile of the bubble wall including CP violation. Here the 
effective potential, which offers the first-order EWPT, should reproduce not 
only Tc and ipc but also various properties of the EWPT such as the latent 
heat and surface tension, which are consistent with lattice results. Since we 
have no available information about the EWPT in the two-Higgs-doublet 
model, we assume some general form of the effective potential at Tc- We 
solve the classical equations of motion for the gauge-Higgs system of our 
model, with the potential replaced by the postulated effective potential, to 
determine the possible profiles of CP violating bubble wall. The equations 
of motion are derived from the lagrangian 

C= --F^^F''^'^ - \b,,B^''' + (D^<l>,)tD^<l>,-Kff(<fi,$2;r), (6.11) 

^ ^ i=l,2 

where 

As long as the EWPT proceeds so that the bubble wall moves keeping the 
shape of the critical bubble with a constant velocity, the static solution con- 
necting the two phases will describe the bubble wall profile. When the bubble 
is spherically symmetric or is macroscopic so that it is regarded as a planar 
object, the system is reduced to an effective one-dimensional one. Here we 
assume that the gauge fields of the solution are pure-gauge type 

tg^^^^ix) = d,U2ix)U,\x), t^B,{x) = d,U,ix)U^\x), (6.12) 

where U2 and Ui are elements of SU(2)l and f/(l)y, respectively. This as- 
sumption would be justified because l+l-dimensional gauge theories contain 
no dynamical degrees of freedom of gauge fields. Further if we find a solu- 
tion based on this assumption, it will have the lowest energy, since it has no 
contribution from the gauge sector to the energy. Then we gauge away all 
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the gauge fields, or, in other words, we fix the gauge in such a way all gauge 
fields vanish. Assuming that U{l)em is not broken anywhere, the classical 
Higgs scalars are parameterized as 







.iOi (x) 



(i = l,2) 



(6.13) 



The static one-dimensional equations of motion are 

2 



- PAZ) 



. dz , 
. d9i{z) 
dz 



eff 



9 / \ dOAz) 9 , , 



dpi 
de2{z) 



0, 
0, 
0, 



(6.14) 
(6.15) 
(6.16) 



dz dz 

where z is the coordinate perpendicular to the planar bubble wall and the 
last equation is the consistency condition for the pure-gauge ansatz, which 
may be viewed as the gauge-fixing condition. We introduce a dimensionless 
and finite-range variable by 



y = -(1 -tanh(az)) , 



(6.17) 



where a has a dimension of mass and its inverse characterizes the width of 
the wall. In terms of this variable, the equations of motion are written as 



d 

4:a'y{l-y) — 
dy 



dpi{y) 

dy 



^a'yil-y)^ 
dy 



1/(1 - y)pl{y) 



V dy , 

de.jyy 

dy 



eff 



dpi 



0, 



OK 



eff 



dOi 



0, 



(6.18) 
(6.19) 
(6.20) 



dy dy 

In order to solve these equations, one must know the explicit form of Vcs- 
Because of the gauge invariance, V^ff is a function of 9i — ^2- From this fact. 
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( 6.19| ) with i = 2 is automatically satisfied as long as p, and 6i satisfy ( |6.19| ) 
with i = 1 and ( |6.20[ ). We determine the effective potential by postulating 
that it is a gauge- invariant polynomial of pi, p2Cos6 and p2sin^ up to the 
fourth order, and that it has two degenerate minima, each corresponding to 
the symmetric and broken phase. Further, for simplicity, we assume that 
the modulus of the Higgs field, Pi{z), has a kink shape and that pi and p2 
have the same order of width. This assumption may amount to that of the 
EWPT proceeding smoothly and accompanying the two Higgs fields which 
take nonzero values at about the same temperature. We expect that the 
effective potential of the model is equipped with these features if the EWPT 
it predicts is first order. 

For the time being, we consider the case in which CP is not explicitly 
violated in Ves, so that all parameters are real and Ves depends on 9 only 
through cos 6'. We require that (|6.18|) has the kink-type solutions in the 
absence of CP violation 

p^{y) = v^{l-y), (6.21) 



where 

^1 = vcos P, V2 = vsin [3. 
Then the effective potential takes the form|[7^ 



Kfr(pi,P2,6') 

{2a? — -mg tan I3)p\ + (2a^ — -mg cot I3)p\ + m\pip2 cos 6 



- {M + 



-2 A cotp + D tan^ p + 



' ' ■ (3 - ' 



fsin/3 cos^/3' 



PiP2(cos( 



+ 



A cot^ P-2DtanP + 



4a2 



(3-- 



vcosP sin P 



PipI{cos9) + Dpi 



A2 4 A3 — A4 2 2 A5 2 2 



^ P1P2 - -JP1P2 cos(2^) 



1 



3 , ^ A2 n \ n 8(3^ 

-Alcot/:/ - — tan /:/ + Astan/:/ . 

2 2 f^smpcosp cos^p 



;(4- 



(6 



P1P2 
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Ai 



cot^ 13 + -A2 tan /3 + A3 cot 13 



;(4 



Here cos^ in the p^-terms is optional. All the parameters should be regarded 
as those including radiative and finite-temperature corrections. Hence even 
for the MSSM, these parameters could be induced in the presence of the 
soft supersymmetry breaking terms. With this potential, the equation for 9 
derived from (|6.19|) in the kink-background is 



P1P2 \ cos 9. 



y) 



• d'9{y) 



+ y{l-y)il-Ay) 



dm 

dy 



d 



[b + c(l - vY - e(l - y)] sm9{y) + -(1 - y)' sm{29{y)), (6.23) 



where the parameters are defined by 



4a^sin (3cos (3 



V 



32a2 
,2 



(Ai cot^ /3 + A2 tan^ /5 + 2A2 



1 



2sin^/3cos2/3 



■^^-^(Ae cot f3 + Xj tan/5), 



d = 



V 

8a 
W 
Aa? 



V 



4a^ sin^ [3 cos^ (3 
V ( B C 



A cos^ j3 + D sin^ (3 



V 



+ 



(6 24) 

4a2 \smP ' cos/5/ ■ 
We refer to the solution of this equation as that with the 'kink ansatz'. When 
9{y) becomes of 0(1), actual solution will no longer satisfy the kink ansatz 
for p{y), ( |6.21| ). Now we show some solutions found numerically for various 



boundary conditions W^. Similar attempts were made by Cline et al.WG 



6.2 Solutions with spontaneous CP violation 

The possible boundary conditions on 9{y) depend on the parameters b, c, d 
and e. Since we concentrate on the potential without explicit CP violation. 
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the boundary condition in the broken phase {y = 0) is either spontaneously 
generated 6q = 6{y = 0) ^ or 6 = mr. The former case is reahzed when 
the parameters satisfy 

d<0 and \b + c-e\<-d, (6.25) 



which corresponds to (|6^ ). Then the boundary value 9q is given by 

cos^o = (6.26) 

Although one might consider that 6i = 6{y = 1) can take any value in the 
symmetric phase since the Higgs scalars vanish there, the finiteness of the 
energy density of the solution requires that 6i = mr. 

We found several solutions, first assuming the power series solution, for 
these two types of boundary conditions. This assumption, which is not es- 
sential, somehow restricts the parameters in the potential. We present some 
of the numerical solutions. 

(i) solution violating CP spontaneously in the broken phase 

When the condition ( |6.25| ) is satisfied, there could be a solution with the 

boundary condition Oq ^ mr and 6i = mn. We found several solutions and 



two of them are depicted in Ref. [Q. Both of them satisfy = 0, while one 
of them has small 6o ^ 0, which may be consistent with the present experi- 
mental bound. The other has 6q = 1. Such a solution may be realized when 
the finite-temperature effects enhance spontaneously generated CP asymme- 
try, which is restored at zero temperature as shown in the MSSM|]7T]. But 



the parameter space admitting such a possibility seems rather restricted!] 
so that the former case with small 6o may be more likely to occur. We show 
the profile for (&, c, e) = (3, 7, 7), 9o = 0.002 and d is determined by ( |6.26| ) in 



Fig. H and the chiral charge flux for the profile at = 0.58 and T = lOOGeV 



in Fig. |10| calculated by the numerical method of Ref. This is almost 
a straight line, because for small 6*0, 9{y) = 6q{1 — y) is an approximate 
solution to the linearized version of ( |6.23|) . As shown in Fig. 10, the magni- 
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Figure 9: The numerical solution of 9{y) for ^(0) = 0.002 and ^^(1) = 0. The 
parameters are (6, c, e) = (3, 7, 7). 
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Figure 10: Contour plot of the chiral charge flux, normalized as 
logj^Q [—Fq/{vwT'^{Ql — Qr))] for the profile shown in Fig. 9. Here we take 
= 0.58 and T = lOOGeV. 
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tude of the chiral charge flux is reduced to about O(10~^) of that in Fig. 
which corresponds to the postulated proflle ( |5.27| ) with maximal CP viola- 
tion. Hence, in this case, even if the forward-scattering is enhanced to give 
an extra factor of O(IO^) so that the baryon asymmetry is the optimal value 
( |5.22| ), the parameter space is rather restricted to generate the present BAU 
only with the charge transport mechanism, 
(ii) solution conserving CP in the broken phase 

In this case the boundary condition are 9q = nn and 6i = mvr. Then there 
are the trivial solution 6{y) = rnr with the kink-type p{y). Besides these, we 
found an interesting solution which violates CP in the intermediate range 
near the bubble wall. If the effective potential admits the condition for the 
spontaneous CP violation ( |6.8| ) to be satisfied for intermediate Vi, as shown 
in the contour plot of Vcs in Ref. ||7^, such a solution exists. Note that 
this condition is weaker than ( |6.25|) , which is required for spontaneous CP 
violation in the broken phase region. Both conditions need A5 < 0. Unless 
the model has tree- level negative A5, it would be difficult to obtain negative 
A5 starting from A5 > at the tree-level, since only the fermions such as 
the gauginos in the MSSM contribute negatively to A5, while contributions 
from the bosons and fermions at finite temperature are positive. We present 
a solution of this type found in Ref. and the chiral charge fiux for it, in 



Fig. |Tl| and Fig. 0, respectively. This type of solution may exist for large 
e, which is the coefficient of p^cos6'-term in the effective potential. As we 
noted earlier, the p^-term arises from the boson loops whose mass at = 
vanishes. In the MSSM, the squark loop might give such a term if its mass is 
very small for (f = 0. Although the maximum of 6{y) is about 0.3, the flux 
is comparable to that in Fig. |^, so that this proflle could generate sufficient 
baryon asymmetry. Since this solution and the trivial solution 6{y) = 
satisfy the same equation and the same boundary conditions, both would 
appear at the EWPT, but with different probabilities. One can determine 
the difference in the nucleation rate by comparing the free energy of the 
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Figure 11: The numerical solution of 6{y) for 6{0) = 6{1) = 0. The parame- 
ters are (6, c, d, e) = (3, 12.2, -2, 12.2). 
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Figure 12: Contour plot of the chiral charge flux, normalized as 
^o§io \~PqI{^wT'^{Ql — <5i?))] for the profile shown in Fig. 11. Here we take 
= 0.58 and T = lOOGeV. 
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critical bubbles. The energy density per unit area of a bubble is given by 

1 




1=1,2 



dz ) \dz J 



dy { ay{l-y) ^ 



i=l,2 
1 



2ay{l - y) 

The difference in the energy density is[[74| 



,dy) ^^'\dy) 
Vcfr(pi,p2,6') \ ■ 



+ V;ff(pi,p2, 



'de, 



/\£ = 8- = -2.056 x 10~^ av^ sin^ (3 cos^ (3. 



(6.27) 



(6.28) 



For the critical bubble of radius Rc, the CP- violating bubble will be nucle- 
ated with probability larger than the trivial one by the factor 



exp 



T, 



c 



(6.29) 



According to the estimation in the massless two-Higgs-doublet model ||39 



the radius of the critical bubble is given by y3Fc/{4:'n-av'^), where Fc is the 
free energy of the critical bubble and is found to be about 145T. Then the 
exponent in ( |6.29|) is 0.89 sin^ /5cos^ /?, irrespective of the wall width. For 
tan/? = 1, the CP-violating bubbles are created about 1.22 times more than 
the trivial ones. 

Now we comment on the baryogenesis in the absence of explicit CP vi- 
olation. Since the effective potential is an even function of 6, —0{y) is also 
a solution to the equations of motion if 6{y) is. The energy densities of the 
both bubbles with 9{y) and —9{y) — we refer to them as 'positive bubble' 
and 'negative bubble', respectively — are equal to each other. Thus their nu- 
cleation rates are also the same, so that the net generated baryon number will 
be zero on the average. This degeneracy of the energy density will be resolved 
once explicit CP violation is taken into account. It was shown that in the 
MSSM, the explicit CP violation in the soft supersymmetry breaking terms 
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induces a mg-type term proportional to sin 6' in the effective potential If 
the difference in the free energy of the two bubbles is AF, the ratio of their 
number will be 

^ = e-^^/^, (6.30) 

and the net baryon number is 

ub (nB\ N+-N_ (nB\ AF 

where {nB/s)^ is the baryon asymmetry generated by each bubble. Comelli 
et al. showed that the explicit CP violation, which does not yield the neutron 
electric dipole moment beyond the experimental bound, splits the free energy 
as large as AF/T ^ 10^^~~^[^. If this is the case, the BAU is explained by 
the electroweak theory as long as {nB/s)o > 10^^. 

6.3 Solutions with explicit CP violation 

As we noted, explicit CP violation is necessary to have nonzero baryon asym- 
metry. We found that even if it is very small, it nonperturbatively yields an 



energy gap between the positive and negative bubbles ||75[|. 

Here we consider only the phase in the m|-term of the type m|(e~*''$|$2 + 
h.c), which is indeed induced in the MSSM. Then the equation for 6 with 
the kink ansatz is now 



r (1 - yy^^ + ^(i - - 4^)^ (6-32) 

= 6sin(5 + e{y)) + [c(l - yf - e(l - y)] sme{y) + ^(1 - yfsm{2e{y)), 

where the parameters are defined by ( |6.24| ). Just as the case with sponta- 
neous CP violation above, the boundary conditions are determined by the 
parameters in the potential. We found two types of numerical solutions. 

One type of solution is that of d{y) = 0{S) in the whole region and 
^(1) = —6. Note that —6{y) is no longer a solution to ( [6.321 ). Hence there 
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is no cancellation, but the generate baryon asymmetry would be the same 
order as that shown in Fig. [TO, if 5 = 0(10^'^). 



The other solution is similar to that in Fig. |TT|, which connects 6{0) ~ S 
and 9{1) = —6 and grows to 9{y) ~ —0.6 at y ~ 0.5. For sufficiently small S 
(~ 10"'^), it has a partner satisfying the same boundary condition but with 
opposite sign in the intermediate region. For larger 6, the partner can no 
longer exist. We calculated their energy density and found that for 6 = 0.0025 
and (6, c, d, e) = (2.98005, 12.178375, -2, 12.2), 

AS = S[e-] - S[e^] = -1.917 X IQ-'^av'^ sin^ (3 cos^ (3, (6.33) 

where 9~^^\y) denotes the solution with lower (higher) energy density. They 
are depicted in Fig. |13], which shows that 6^{y) deviates from the solutions 
with (5 = by 0(0. 1) near the bubble wall {y ~ 0.5) in spite of the small S. 
The ratio of the nucleation rates of the bubbles is 

- = exp ^ j = 8.05, (6.34) 

where we put tan/3 = 1. Then the net generated baryon asymmetry will be 

^"-'^ +('^) -J^^fHS.) (6.35) 



s \ s J + N+ + N_ \ s J ^ N+ + N_ \ s J + + 
where we used (ub/s)- ~ — (nB/s) + .[] Since the chiral charge flux for this 



profiles are expected to be the same order as that in Fig. |T^, we will obtain 
sufficient baryon asymmetry. 

All the solutions presented in this section are based on the kink ansatz. 
In fact, the solutions other than the trivial ones 6{y) = mr are not exact 
solutions to the full equations of motion. We also get solutions without this 
ansatz, which would have lower energy. In fact we found some numerical 



solutions, which have 9{y) of 0(1) so that Pi{y) are no longer kink shape |]77| 



'^When the discrepancy between 0'^{y) and —9 (y) is large, this approximation is in- 
vahd. 
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Figure 13: The numerical solutions 9 (y) (the dashed curve) and 9 (y) (the 
dotted curve) for 5 = 0.0025 and {b,c,d,e) = (2.98005, 12.178375, -2, 12.2), 
which satisfy the boundary conditions, 9^ ~ —0.3115 and 9f = —5. 
The sohd curves are the degenerate solutions for 5 = and (6, c, d, e) = 
(3, 12.2, -2, 12.2) with = = 0. 
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Such solutions may exist for a broader parameter region. We have concen- 
trated on static solutions. If the EWPT accompanies exploding bubbles, we 
would need to solve the time-dependent equations, and if the viscosity of the 
medium cannot be ignored, we would have to solve the equations with fric- 
tions or random forces and then the wall profile might not give a monotonical 
p like those obtained here. These extensions of the equations of motion are 
to be studied in the future. 

7 Discussion 

In this paper, we briefly reviewed the scenario of the electroweak baryogen- 
esis. For it to be successful, some extensions of the MSM would be needed 
to have CP violation in the Higgs sector. As one such model, we analyzed a 
two-Higgs-doublet model, including the MSSM, and tried to relate CP viola- 
tion to the baryon asymmetry generated at the EWPT. The main attraction 
of the electroweak baryogenesis is that it relies exclusively on physics which 
can be tested by present or near future experiments. In order to bridge the 
gap between detectable microphysics and the baryon asymmetry observed in 
the universe, we still have to know various aspects of the electroweak theory, 
in particular the EWPT and CP violation in the model. 

In principle, one can predict how much the baryon asymmetry is gen- 
erated at the EWPT, once the renormahzed lagrangian of the electroweak 
theory is given. Prom the lagrangian, one can construct the effective potential 
at finite temperature near the phase transition. One may use higher-order 
perturbation or nonperturbative methods such as the lattice simulation. As 
we saw in § 4, the studies have been limited to those of the MSM. We hope 
these efforts to be extended to the extended versions of the MSM. The ef- 
fective potential gives information about the EWPT, such as its order, the 
transition temperature, and the latent heat and surface tension if it is first 
order. If the EWPT is not first order, the baryon number would be washed 
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out except for the portion proportional to the primordial B — L. Even if 
the EWPT is first order, the unsuppressed sphaleron process would erase the 
baryon asymmetry. The Higgs mass is bounded from above, requiring that 
the EWPT be first order or that the sphaleron processes decouple after the 
EWPT. The upper bound seems inconsistent with the present lower bound 
of the Higgs scalar in the MSM. More detailed studies of the EWPT in the 
extended models will provide new constraints on the lightest Higgs boson. 
When the EWPT is first order, one can describe how the transition proceeds. 
If the phase transition accompanies nucleation and successive growth of the 
bubbles of the broken phase in the symmetric phase, the mechanisms pre- 
sented in § 5 will yield the baryon asymmetry. The estimate of the generated 
baryon number depends on the velocity and thickness of the bubble wall and 
CP violation around it, both of which should be dynamically determined. 

When we analyzed the equations of motion for the Higgs fields in § 6, we 
did not completely fix the parameters in the potential. In practice, these are 
fixed once one chooses the model so that the type of profile is realized would 
also be determined. Under the assumptions made there, A5 < might be 
necessary to have large CP violation, while it may be a weaker condition for 
spontaneous CP violation in the broken phase. Then a very small explicit 
CP violation is needed to have nonzero net baryon asymmetry. In any case, if 
more detailed features of the EWPT in the two-Higgs-doublet model become 
available, one could calculate the generated baryon asymmetry, which in turn 
constrains some parameters in the model to explain the present BAU only 
with the help of the electroweak theory. 
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A Summary of the Big Bang Cosmology 

In this appendix, we summarize the basics of the standard big bang cosmol- 
ogy and estimate various time scales at temperature near the electroweak 
phase transition. The most general form of homogeneous and isotropic space 
is the Friedmann-Robertson- Walker (FRW) metric, given by 



ds' = dt'-R\t) 



1 — fcr^ 



(A.l) 



where R{t) is the scale factor in the comoving coordinate, k characterizes the 
topology of the space and k = 1,0,-1 correspond to closed, fiat and open 
space, respectively. The Einstein equation leads to 

Uj +^-3 = ^^-2) 

|-f - + (A.3) 

where p is the energy density, p is the isotropic pressure, and A is the cos- 
mological constant, p and p are related by the equation of state; p = ■yp, 
where 7 = 1/3 for a radiation-dominated universe and 7^1 for a matter- 
dominated universe. Eliminating A from (|A.2|) and (|A.3D, we have 



{R^p} + m^Rp = 0, (A.4) 
from which we obtain by use of the equation of state 

p^3{7+i) ^ constant. (A.5) 



Because of 



8 30 ' 
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the energy density in a radiation-dominated universe is given by 

p{T) = ^g^T\ (A.7) 
where represents the effective degrees of freedom at temperature T, 

9* = T.9b + IT.9f. (A.8) 

B ^ F 

Here gB{F) counts the spin and internal degrees of freedom of the bosons 
(fermions). For the standard model with Nf generations and m Higgs dou- 
blets, 

7 

c/, = 24 + 4m + - X 30Nf, (A.9) 
8 

so that g^ = 106.75 for the MSM. In a radiation-dominated universe, the 
Hubble parameter is approximately given, from ( |A.2| ) with A = 0, by 



i/^j^p^ 1.667^ A (A.IO) 
V 3 mpi 



where mpi = \/Gn = 1-22 x lO^^GeV is the Planck mass. 

The entropy density at present is related to the photon density as follows. 
As long as the local equilibrium is maintained, the law of thermodynamics 
implies 

d{p,y) = TidSi - PidV + pidNi, (A. 11) 

where the index i denotes the particle species and Si is the entropy in a 
comoving volume V . This equation, together with the Gibbs-Duhem relation 

SidTi - Vdpi + Nidfii = 0, (A. 12) 

lead to, up to a constant, 

^ p^ + p^ + pj^_ 
Neglecting the chemical potentials (/ij ^ Tj), the entropy density is given by 

i i * 
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Since the entropy is dominated by relativistic particles 



45 

where 



9.s^^ (A. 15) 



We can ignore the difference between Tj and T at high temperatures as 
shown below, while the difference is significant for T < IMeV. At such 
low temperatures, v-v annihilation decouples and the entropy in e^ pairs is 
transferred only to the photons. The entropy conservation implies 

= + = (A.17) 



Hence today we have 

(7,5 = 2 + ^x^x2x3 = 3.91, (A.18) 

o ii 

which yields 

s = 7.04-n^ (A. 19) 

with being given by ( |A.21| ). 

Now we shall estimate time scales of various interactions near Tc — 
lOOGeV of the EWPT. Given a cross section a of some interaction, the mean 
free path A is estimated as 

A-cT = -, (A.20) 

n 

where n is the density of the particles which participate the interaction. At 
temperature T, n of a massless particle whose degree of freedom is g is 



d-k I I ^-^gT 



^3 



^ = ^/ (2^)3 eN/^ + 1 = 1 3^03) ^3 ^^-''^ 

4 7r2 ^ ' 
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where ({3) = 1.2020569 ■ • ■. The mean free time t of a particle of mass m 
and energy E is given by 

t = -= , ^ (A.22) 

For E :> m, t c::^ X. Since the cross section of the interaction with the fine 
structure constant a at the center-of-mass energy -^i is o" ~ a'^/s, the mean 
free path at T is 

10 1 
^ ^ ^ ■ - ^ TTT^' (A-23) 

where we have used g^, ^ 100 and s ~ T^. 
If we take T = lOOGeV, 

As ~ — - — ~ O.lGeV"^ for strong interactions, 

Xew — — 5 ~ IGeV"^ for electroweak interactions, (\ 24) 

lO'^av^ ^ • ^ 

Ay ~ I — — ) Xew for Yukawa interactions, 

J 

where we have used as{mz) = 0.117±0.005 and aw = C(qed/ sin^ 9w — 1/30. 
At this temperature, the time scale of the universe expansion is, from ( [A.10| ), 

H-\T) ~ lO^^GeV-^ (A.25) 

The time scale of the sphaleron process is 

tsph ^ (rsph/ri)-' ~ 10^GeV-\ (A.26) 

The thickness and velocity of the bubble nucleated at the EWPT are 

1 ~ 40 

~ ^ - 0.01 ~ 0.4GeV-^ (A.27) 

and 

Vy, ~ 0.1 ~ 0.9, (A.28) 
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respectively [^]. Although and correlated, we adopt as the scale of 

the EWPT, 

t^^ii = hLr^ 0.01 ~ AGeV-\ (A.29) 
From these we see the following: 

1. All the particles are in kinetic equilibrium at the same temperature, 
because of H^^ » tEw, far from the bubble wall. 

2. The Yukawa interactions of the light fermions {nif < 0.1 GeV) are out 
of chemical equilibrium. 

3. Some of the flavor- changing interactions are out of chemical equilibrium 
because of small KM matrix elements. 

4. Since for the leptons Ay > \ew ^ ^w, the leptons propagate almost 
freely before and after the scattering off the bubble wall. 

5. Because of twaii ^ isph, the sphaleron process is out of chemical equi- 
librium near the bubble wall. 



B Derivation of the Baryon-Number-Changing 
Rate 

Suppose that there are states with discrete labels, i. Let P{i; t) be the 
probability to find the system in the state i at time t, and Ti^j be the 
transition probability from state i to j per unit time. Then the following 
master equation holds: 

Pit; t + At) = -Y. Pit; t)r,_, At + J] P(j; t)T,^,At + Pit; t). (B.l) 

For a steady state, detailed balance is maintained, that is, P(i; t) is indepen- 
dent of t; P(i; t) = Pcqii) for any t. Then ( p.l|) is reduced to the detailed 
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balance equation; 

This implies for the case of baryon number changing processes 

oo 

+ rB 

n=l 

oo 

= Y.\P-^(B + n)TB+n-,B + P.^{B-n)VB-n^B\. (B.3) 



n=l 



where Pq^{B) oc e~^^/'^ with Fb being the free energy of the state with baryon 
number B. The transition rates with Ai? = ±1 are given by r+ = T b^b+i 
and r_ = Vb^b-i and the others are approximated by TB^B+n — T" and 
^B^B-n - . Since Fb+„ = Fb + rifiB, (|B.3|) yields 



OO OO 

E [r+ + ~ E [(e'^^/^r_)" + (e'^^/^r+)"] . (b.4) 

n=l ?i=l 

If F-t ^ 1, only the contribution from = 1 in the sum is dominant so that 

r+ + r_ ~ e"^'s/Tp_ _^ e''s/^r+, which leads to 

£± ^ e-^^/'^. (B.5) 



By definition, = F^ — r_ if we use r± as the rate per unit volume and unit 
time. That is approximately given by the sphaleron rate, r+ ~ r_ ~ Tsph- 
Hence we have 

UB = r_ (^^ - ^ r3pUe---/^ - 1) ^ (B.6) 
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